Effective merging dynamics of two and three fluid vortices: 
Application to two-dimensional decaying turbulence 
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We present a kinetic tlieory of two-dimensional decaying turbulence in tiie context of two-body 
and three-body vortex merging processes. By introducing the equations of motion for two or three 
vortices in the effective noise due to all the other vortices, we demonstrate analytically that a 
two-body mechanism becomes inefficient at low vortex density n <^ 1. When the more efficient 
three-body vortex mergings are considered (involving vortices of different signs), we show that 
n ~ , with 5 = f. We generalize this argument to three-dimensional geostrophic turbulence, 
finding ^ = 5/4, in excellent agreement with direct Navier-Stokes simulations [J. C. Mc Williams et 
al, J. Fluid Mech. 401, f (i999)]. 

PACS numbers: 02.50.-r, 47.10.-g, 47.27.-i 



I. INTRODUCTION AND BACKGROUND 

Two-dimensional freely decaying turbulence (i.e. when 
an external forcing has been switched off at the initial 
time t = 0) exhibits fascinating properties starting 
from an initial incoherent turbulent background jTobust 
coherent structures (the vortices) soon emerge Dif- 
ferent types of structures can be observed: monopolar 
vortices rotating in either direction lZH9| , quasi ballistic 
dipoles (a vortex/antivortex pair) [9|, llOl ] and even rare 
tripoles (a pair of two like-sign vortices separated by an 
antivortex) [l^ . These vortices interact due to their mu- 
tual advection which amounts to a long-range interaction. 
The vortices are characterized by their radius a{t) and by 
their core (or peak) vorticity uj{t) while their precise vor- 
ticity profile is not particularly relevant in the process of 
2D decaying turbulence. When two like-sign vortices ap- 
proach at a critical distance of the order of their radius 
a{t), a complex merging process occurs pT| and results 
in the formation of a larger vortex. As a consequence 
of the merging processes, the vortex density n{t) decays 
with time while their typical size a{t) increases. 

Numerical simulations of the Navier-Stokes (NS) equa- 
tion [Si, 12 18j and of the point vortex model presented 
below [ij, [13] are consistent with experiments [2l| - 
[2^ in finding a power-law decay of the vortex density 
n ^ t~^, with ^ in the range 0.6— 1.0, although it is most 
often found around £, ~ 0.70 — 0.75. For an inviscid fiow, 
the dimensional Batchelor argument [2^ which assumes 
that the energy per unit area, E ^ i^/L'^) J rf^x ^ 
uj'^na^, is the sole invariant predicts n ~ l/{Et^). This 
implies ^ = 2, in clear disagreement with numerical and 
experimental data. However, in [l^ . and in agreement 
with subsequent theoretical and experimental work, the 
average peak vorticity uj was also found to remain almost 
constant during the dynamics (see Appendix IX| . 

The theoretical determination of the exponent ^ has 
been a subject of active research and many theories. 



based on different physical arguments - and leading in 
general to different results - have been proposed over 
the years [H, [13, [26l - [30j . Recently, [3l| reviewed several 
theories and showed that some of them are inconsistent 
(except that of [13] and [13l; [13 "^^^ not reviewed, since 
not yet published at the time) with the fact that the 

,r7l/2„2 , 



vortex typical velocity v ^ cjn^^'^a'^ ^ V E is constant 

The vortex density satisfies a general equation of the 
form 



dn/dt = ^n/T„ 



(1) 



where the merging time T,„(i) is, by definition, the typi- 
cal time needed for two vortices to actually merge. Since 
we expect a power-law behavior for n{t) ~ t~^, we con- 
clude that T„i ^ t. Note that checking that t,„ ~ t in 
NS simulations or experiments is a good test to know 
whether the system has entered the asymptotic scaling 
regime [13] • In units of the collision time. 



R/v ^ {ojna^) 



(2) 



which is the typical time needed for a vortex of veloc- 
ity V ~ Lon^^'^a'^ [201 to travel the average inter vortex 
distance R ~ n~^'^, we obtain the most general form 



R 

T I — 

a 



(3) 



where the exponent a cannot be determined on purely 
dimensional grounds since na^ ot R/a are dimensionless. 
The area fraction covered by the vortices na^ is actually 
the only available dimensionless quantity [13]. Now im- 
posing energy and peak vorticity conservation, Eq. ^ 
leads to 
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(4) 
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SO that determining ^ amounts to determining a, or the 
merging time r,„. Note that Batchelor's argument [1^ as- 
sumes that Tm ^ T ^ t, i.e. a = 0, after which the energy 



2 



conservation and Eq. ^ leads to ^ = 2, and a ^ R ^ t, 
in clear disagreement with numerical and experimental 
data. Hence, the necessity of a more general scaling the- 
ory arises from the fact that the two length scales a (typ- 
ical vortex radius) and R (typical inter-vortex distance) 
do not scale in the same way with time. In other words, 
the mean free time r and the merging time arc not 
equivalent quantities, especially at low vortex density. 

We now consider the effective vortex model introduced 
in p^ . The exact Hamiltonian for point-like vortices 
can be derived from the inviscid NS equation [32| and is 
reminiscent of two-dimensional electrostatics 



U = -^7,7jln(|ri 



(5) 



where 7^ ^ ±wa^ are the vortex circulations. The crucial 
difference with electrostatics resides in the fact that the 
vortex coordinates {xi,yi) are the conjugate variables, so 
that the resulting equations of motion are 



dxi _i dH _ Hi - Hj 



dyi 
dt 



dxi 



(6) 
(7) 



Assuming that all vortices have the same conserved peak 
vorticity iw, these equations are supplemented with the 
following merging rules consistent with NS simulations 
113I Il9| and experiments: when two like-sign vortices 
of radius ai and 02 reach a distance Vc ^ ai 4-02, they 
merge keeping the energy E ~ nuj'^a'^ constant, so that 
the resulting vortex has a radius a satisfying, = af + 
a\. The average radius then scales as a ~ {nitP' j E)~^ I ~ 
t^l'^ R, where R ~ ^ t^/^ is the typical distance 

between vortices. Hence, for large time for which a <^ R, 
we expect that the application of the point-like vortex 
model to two-dimensional decaying turbulence becomes 
more and more justified. 

The different scaling relations obtained in this work 
rely on the assumption that the vortices are character- 
ized by a single length scale a. Hence, our results will 
be valid provided the vortex radius distribution p(r, t) 
is integrable or marginally integrable for r 0, i.e. 
p{r,t) <^ , with 6 < 1. Otherwise, when the vor- 
tex size distribution is polydisperse (for S > 1), one must 
introduce a cut-off r< for small radii in order to ensure 
that p{r, t) is normalizable. Then, r< and a will scale 
differently with time, thereby introducing another rele- 
vant length scale apart from a and R. In the case 5 = 1, 
r< and a scale identically, up to logarithmic corrections 
in the time and the present theory still applies. Start- 
ing from an incoherent initial background with an initial 
energy spectrum peaked at a certain scale fco = l/oo, ex- 
periments and NS simulations show that the radius dis- 
tribution remains integrable and is very often bell-shaped 
(in the NS of [l^, it is however claimed that 5 = 1) and 



scale-invariant, hence satisfying 



(8) 



where g is the scaling distribution and the prefactor en- 
sures the normalization of p. Similarly, the peak vorticity 
distribution is found numerically and experimentally to 
be bell-shaped, so that only a single vorticity scale w 
must be introduced. 

From now on, we place ourselves in the framework 
of the point vortex model described above, for which 
we present analytical and numerical results. Hopefully, 
and as claimed in the two preceding paragraphs, this 
model is faithful enough to the original problem to ap- 
ply our results to decaying turbulence in actual fluids. 
Our first aim in Section [H] is to evaluate analytically and 
numerically the time - hence the exponents a and 
^ - necessary to observe a merging in the context of a 
purely two-body mechanism. In Section IlIIl we will then 
address the relevance, at least at small vortex density 
na^ ^ 1, of three-body mergings (typically, a ballistic 
dipolc hitting a single vortex), which will be shown to ul- 
timately dominate. The numerical results of this second 
study will be supplemented with a fully consistent ana- 
lytic mean-field three-body kinetic theory showing that 
the asymptotic density decay exponent should be ^ = 1. 
We finally apply a similar analytic argument to three- 
dimensional geostrophic decaying turbulence, finding an 
exponent ^ = 5/4 in excellent agreement with direct NS 
simulations [ssj . The present work hence precisely as- 
sesses the crucial differences between a vortex dynam- 
ics dominated by two-body or three-body (a dipole and 
an isolated vortex) merging processes, demonstrating the 
leading role of the later at low vortex density (the impor- 
tance of three-body interactions in sufficiently dilute 2D 
turbulence was mentioned in [sij and incorporated in ki- 
netic models of decaying 2D turbulence in [3 113] , albeit 
in a different manner). 



II. THEORY OF TWO-BODY MERGING 
PROCESSES 

We first consider the effective dynamics of two inter- 
acting like-sign vortices, i = 1,2, which are also ad- 
vected by the other vortices, assumed to be uniformly 
distributed in space but remaining at a distance larger 
than R from both tagged vortices (two-body merging 
mechanism). Since we are ultimately interested in the 
regime r = |ri — r2| ^ a <^ R, one can expand Eqs. (|6l7p 
in powers of the relative distance between both vortices, 
r. After expressing the time in unit of r and distances 
in unit of i? ^ 77,"^/^, we find that the relative distance 
between both vortices r = ri — r2 obeys 

dr 

r[cos{2ip)ria + sm{2(p)rip], (9) 



dt 

dip 



1 



sin(2(/))?7Q -I- cos(2(^)77,3. 



(10) 
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where we have used cyhndrical coordinates r = (x, y) = 
(r, tp) and the fictitious particle at r is now confined to 
a disk of radius R ~ 1. Furthermore, r]a and rjp repre- 
sent the effective noise due to all other vortices, and are 
explicitly given by 



'7a = ^3 



2xojyoj 



VP 



(11) 



where rg = {xo,yo) is the position of the center of mass 
of the two vortices. In the absence of the other vor- 
tices {r]a = rip = 0), the two tagged vortices would ro- 
tate around each other remaining at a constant distance 
r. Moreover, even in the presence of the other vortices, 
dr/dt — J- as r — 7> 0, since if both vortices were at the 
same position, they would experience the very same ad- 
vection from the other vortices. The noises can be shown 
to satisfy ^ 



ivlit)) = ivlit)) 



(77„(t)r7^(0) -0, (12) 



where ~ C'(l) in onr dimcnsionlcss units. Moreover, 
these noises only vary notably on the time scale over 
which the other vortices travel a distance of order R. 
Their correlation time is thus of order r ^ 0{1) [1^. The 
fact that this correlation time is non-zero (as it would be 
for a standard pure white noise) is crucial: when the 
two vortices arc very close to each other, their mutual 
rotation period will become much smaller than r, which 
will result in the effective averaging out of the external 
noises r]a^p. Finally, it appears reasonable to modelize 7?a 
and rjfj by effective Ornstein-Uhlenbeck processes [lO, HI] 
with variance and correlation time both equal to 1: 



dt 



77 



7 = q;,/3. 



(13) 



where Wa and wp are two independent (5-correlated Gaus- 
sian white noises, implying {'r]j{t)VY(^')) ~ '<^7,7'- 
The Fokker-Planck equation associated to the system 
of Eqs. (|9ll0ll3p . which describes the evolution of the 
probability distribution function (pdf) of the variables 
{r,ip^ria,rip), can be easily written (sec Appendix IB|). 
Unfortunately, it seems unlikely that its general time- 
dependent solution can be analytically obtained. How- 
ever, it is straightforward to check that P(r, ip, rjai Vp) 
r e~'^o,/'^-Vfi/'^ jg stationary solution. Hence, the relative 
coordinate r — (x, y) has a uniform asymptotic distribu- 
tion. 

For r <C 1, we can still obtain some useful information 
on P{r,t), since the stochastic term in Eq. pH)) can be 
neglected with respect to in this regime. Dividing 
Eq. ^ by Eq. ([TU]), and integrating over (p, we obtain 



1 



1 



cos{2(p')r]a ((/?') + sm{2(p')7]i3 {ip')] dip' . 

(14) 

Hence, for small r, the variable r ^ has approximately a 
Gaussian distribution (being the linear sum of Gaussian 



variables r]a and ry^) of variance 



d^i d(p2 cos[2(</J2-(^i)]e-l*^^=)-*('^^)l. (15) 



Since t'{(p) « r'^if) ^ 1, we can write t{(p2) — tifi) ^ 
{ip2 — ipi)t'{ipi). Finally, after computing the integral in 
Eq. ([15]) in the limit of large <j>, 

d^2 cos[2(^2 - ,,0]c-(— = 



(16) 



we find, using Eq. (jl4p . and for large ip and time t 



1 



1 



j-^it) r2(0) 



t 

t'{ipi)dipi 



(17) 



Remembering that the asymptotic pdf is P{r) ^ r and 
that the random variable is approximately Gaussian 
for small r, the above argument suggests that P(r, t) can 
be written as 



P{r,t) = rf[r/roit)l ro{t)^t-^/\ f{x) 



(18) 



where the typical scale ro{t) ^ t ^l'^ has been deduced 
from Eq. ^7)) . 




r/r^Ct) 



FIG. 1. (Color online) Insert: we plot P{r,t)/r for 7 dif- 
ferent times (in geometric progression) up to f ~ 10^, av- 
eraged over trajectories starting from r(0) ~ 1. The main 
plot shows the corresponding data collapse, fully consistent 
with Eq. (|18|l . along with a fit to the predicted Gaussian 
functional form (in the variable x~'^) for the scaling function, 

f{x) = e-(^"'-^(^')' {xo ^ 2.1; dashed line). 

This scaling behavior is illustrated in Fig. [TJ where 
Eqs. (|9ll0ll3p have been solved numerically using a sec- 
ond order stochastic integration scheme. We ran ~ 10^ 
trajectories up to < '--^ 10^, starting from random initial 
values for r(0) '--^ 1 (i? = 1 in our units) and (f)(0) (uni- 
formly distributed in [0, 27r]). The distribution of r{t)^ 
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P{r,t), is plotted for seven different times, illustrating 
the perfect scaling collapse predicted by Eq. p^ . 

Using Eq. (fT8|) . we can obtain the behavior of the mo- 
ments of r{t) for large time. We find that for z > —2, 
mz{t) = {r^[t)) converges to a constant, while m-2{t) ~ 
\a.t, and 



(19) 



In addition, Eq. dTO]) implies {ip{t)) ^ t{r-^{t)) ^ tint. 
These different predicted behaviors are in excellent agree- 
ment with numerical simulations of Eqs. (|9ll0ll3p pre- 
sented in Fig. O 




FIG. 2. (Color online) From bottom to top, we plot m_4(t) = 
{r~^{t)), exp{{r~^{t))) (which should have a power-law be- 
havior, since m_2(i) ~ In^), and {ifi{t)). The dashed lines 
are fits to the predicted functional forms (see text): aot°'^ 
for exp[m_2(t)] and m_4(t) (in this case, ai — 0.49(1)), and 
aotln{t + Oi) for {(p{t)). The Insert illustrates the near con- 
stancy of m-(i{t)/m^4{t), which is also fully consistent with 
ro{t) ~ f"-^/" and the scaling form of P{r,t) in Eq. (fTSl) . 



Let us now address the problem of determining the 
merging time Tm{rc) for two vortices of radius Tc start- 
ing at a distance r(0) ~ 1, both lengths being expressed 
in unit of R. By definition, Tm(rc) is the average time 
necessary for the distance r between the two vortices to 
reach the value Tc for the first time. Alternatively, defin- 
ing rmin(i) as the minimum r(t) reached up to time t, 
one should have 



(^min ('^m (^c) )) ^ 



(20) 



and one can obtain an estimate of rm(rc) by inverting 
this relation. Moreover, our previous heuristic argu- 
ment resulting in Eq. pT)) actually shows that the av- 
erage escape time Te(rc) to go from a region for which 
r(0) e [(1 - e)rc, rd (with < £ < 1) to another for 
which r{t) G [1 — £, 1] is given by 



-ire) 



(21) 



A (Markovian) detailed balance argument then implies 
that the typical time to go from r(0) S [1 — e, 1] to r{t) G 



[(1 — e)rc, Tc] - a time we assimilate to Tm(rc) - satisfies 
Prob(r e [1 - e, 1]) Prob(r e [(1 - e)rc, r^) 



Tm{rc) 

Using the stationary pdf, P{r) 
for Te(rc), we obtain 



r, and the above result 



t 



-1/6 



(23) 



We have numerically estimated {rm\n{t)) and the 
first-passage times T„i(rc) and Te(rc), by integrating 
Eqs. (|9ll0ll3p using a second order stochastic integra- 
tion scheme (sec Fig. being initially fixed, the 
merging time T,„(rc) is determined by averaging over 
^ 10^ (large Tc) to ^ 10^ (smallest Tc) realizations, the 
first time for which the inter- vortex distance r[t) < Tc 
(merging criterion), starting from an initial r(0) ~ 1. 
Conversely, starting from an initial r(0) = Vf. < fmax, the 
escape time Te{rc) is estimated as the first time for which 
r{t) > Tmax ^ 1 (we havc set rmax = 0.8 in unit of R). 
The results presented in Fig. [3] are in good agreement 
with our analytic estimates of Eqs. (|21l23p . 
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FIG. 3. (Color online) For the two-vortex effective model 
of Eqs. (|9I10I13|I . we plot Tm,(r^ vs r^^ (top full line) and 
the time t vs ('■~f^(t)) (middle) along with power-law fits 
with respective exponent a « 3.1 and a w 3.0 (dashed fines). 
We also plot the escape time Te(rc) (bottom) and a fit using 
a — 2. These exponents are in good agreement with the 
predictions of Eqs. P7l21l23p . For the three- vortex model 
(Insert), we plot Tm.{rc) and power-law fits with exponent a « 
1.8 for "large" Tc (high density) and the predicted asymptotic 
exponent a = 1, for rc — >■ (low density). 

Reintroducing time and spatial units, we conclude that 
a purely two-body mechanism would lead to a = 3 (see 
Eqs. (|3I23[) ) and ^ = 1/2. Even if we assimilate the ef- 
fective cut-off ro(t) ^ R(t/T)~^/^ to the vortex radius 
a(t), we end up with a = 2 and ^ = 2/3. We will show 
below that a three-body merging mechanism leads to a 
faster density decay, and hence should dominate the dy- 
namics at low vortex density. Actually, the inefficiency 
of the two-body merging processes has a simple physi- 
cal interpretation suggested by our analytical results of 



5 



Eqs. (|14I15I17P : when two like-sign vortices are at a dis- 
tance r ~ a ^ i?, they rotate around each other at such 
a high velocity {d(p/dt ^ r~^) that the effective noise 
due to vortices at a distance greater than R (which has a 
correlation time of order r ~ 1, in our units) is averaged 
out by the fast rotation. In other words, the fast rota- 
tion makes the effective noise felt by both tagged vortices 
so weak that the two-body system becomes almost intc- 
grable (in the absence of noise, the distance between both 
vortices would remain strictly constant). This points to 
the necessity of the presence of at least another vortex 
at a distance of order a, in order to significantly perturb 
the nearly integrable two-body system. 



III. 



THE ROLE OF TREE-BODY MERGING 
PROCESSES 



By generalizing our model of Eqs. (|9I10I13I) to the case 
of 3 vortices (2 like-sign vortices and 1 of opposite sign), 
we define again r„i(rc) as the time for which the distance 
between both like-sign vortices reaches the value rc for 
the first time. We numerically obtain that Tmirc) ~ r^'^ 
(see below for a theoretical argument), leading to a = 1 
and = 1- The data are presented in the Insert of Fig.|31 
confirming the results obtained in [2^, although for a 
wider range of values of rc- However, we find that the 
high density regime, corresponding to the largest values 
of Tc, is best fitted by a « 1.8 (see Insert of Fig. [3]). This 
value of a is associated to ^ « 0.71, in good agreement 
with (high density) simulations and experiments. More- 
over, we find that when the two like-sign vortices are at 
a distance rc <C 1, the third (opposite-sign) vortex is at 
a distance proportional to rc itself, justifying that this 
third vortex plays a crucial role in the collision process. 
Hence, we have shown that considering three-body merg- 
ing processes speeds dramatically the dynamics, at least 
at sufficiently low densities na^ <^ 1 . This leads to a true 
scaling exponent ^ = 1 (associated with a = 1). This 
regime is, however, difficult to reach in NS simulations 
and experiments, and we claim that only a pseudo-scaling 
is observed at the densities achieved, na^ ^ 0.1. This 
corresponds to a transient regime with a pseudo-scaling 
exponent ^ ~ 0.71 (associated with a ~ 1.8) measured 
on less that one decade. This pseudo-scaling exponent 
is expected to slowly increase with time, asymptotically 
reaching the true scaling exponent ^ = 1. A strong evi- 
dence that true scaling is not reached in the experiment 
of [23} is that the measured merging time was found to 
scale as Tm ^ t0.57±o.i2^ instead of the necessary r™ ^ t. 
In addition, an indication that the exponent ^ measured 
in direct numerical simulations slowly increases with time 
is given in [l3l . In order to reach the true scaling regime, 
one would need to run longer simulations with more vor- 
tices at the start which is difficult. Alternatively, a dy- 
namical renormalization group analysis [20j allows one 
to reach very low densities more easily. In this case, it 
is found that ^ = 1 (measured on four decades) in agree- 



ment with the present three-body kinetic theory. 

Note that the local energy (see Eq. ([5])) involving three 
vortices of respective circulation 7, 7, —7, is ^ 
7^ In(ri3r23/ri2). It can remain finite and even constant 
even if ri2 <^ 1, provided the third vortex of opposite sign 
is in the vicinity of the two merging vortices. In the case 
of a binary merging, the local energy ^ —7^ lii(''i2) 
diverges when ri2 <SC 1 and the excess energy must be ab- 
sorbed by the far away vortices, a very slow process, as 
shown in the present work. A similar argument explains 
why a three-vortex merging process involving three like- 
sign vortices is as inefficient as a two-body process, as 
observed numerically in [20| . Finally, we have confirmed 
numerically (simulation not shown) that adding a fourth 
vortex in our effective model docs not change the scaling 
of the Insert of Fig. [31 provided that all four vortices do 
not have the same sign. 

At low density na"^ ^ 1, the proposed physical picture 
is that ballistic dipoles collide with isolated vortices, the 
later ultimately merging with the like-sign vortex of the 
dipole. We briefly recall the kinetic theory presented 
by two of us in |20[ and generalize it to any dimension 
d, and in particular to geostrophic decaying turbulence 
of a three-dimensional stratified fiuid [33|, where planar 
vortices are replaced by almost spherical "blobs". The 
vortex density satisfies an equation of the form dn/dt ~ 
—n/Tfn and, as already mentioned, the density power-law 
decay n ~ t~^ implies r,„ ~ t. Noting that the dipoles 
move almost ballistically with velocity Vdip ~ 7/a ^ 
and using a classical cross section argument, we obtain 



l/(n* 



Vdip), 



(24) 



where d is the dimension of space. Using a mean- 
field approximation, dipoles of typical size a are formed 
with a density ridip ^ nxna'^. Moreover, the energy 
E ~ iujay^na'^ and the (potential) vorticity peak lo are 
conserved, so that Eq. ([M)) leads to 



E 



xn <i+2 . 



(25) 



where we have kept track of the dimensional constants uj 
and E. Equivalently, this result can be rewritten in the 
same form as in Eq. (|3]) 



(26) 



(na^) 



6-d 
2 \ d+2 



Finally, expressing the asymptotic scaling condition, 
T^^t, Eqs. (|25I26P lead to 



d+2 6-d 



(27) 



In d = 2, we obtain ^ = 1 (associated to a = 1), which 
agrees with the numerical results of Fig. |31 Again, this 
value is slightly higher than the exponents measured in 
experiments and most numerical simulations. However, 
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we have already argued that the smallest area fraction 
na^ reached remains quite large, except in the RG proce- 
dure of [i^l, which indeed found f w 1. In the experiment 
of (23j , the measured merging time was found to scale as 
Tm ^ iO-5''±o.i2^ instead of the necessary Tm ^ t, showing 
that the scaling regime was not yet attained. However, 
we observe that the theoretical relation of Eq. ([25)1 is 
remarkably well verified in this experiment: 



t' 



0.57±0.12 



E 



■xn 



0.55±0.14 



(28) 



This suggests that the result of Eq. (P5|) could be valid 
in a wider range of time, even before reaching the true 
asymptotic scaling regime. This phenomenon is actually 
quite common in the context of critical phenomena, when 
dealing with finite-size samples and simulation time: the 
scaling relations between critical physical quantities (like 
the susceptibility or the correlation length/time in mag- 
netic systems) are often found to be better obeyed nu- 
merically than between these quantities and the critical 
parameter itself {e.g. the distance to the critical temper- 
ature |T-Tc|). 

In d = 3, we obtain ^ = 5/4 (associated to a = 3/5), in 
good agreement with the value ^ ~ 1.25 ± 0.10 measured 
in NS simulations of gcostrophic decaying turbulence in 
a stratified fluid [11]. We argue that due to the fastest 
decay in d = 3, the low density regime where the three- 
body mechanism sets in appears faster. In addition, the 
increased dimensionality generally plays in favor of a bet- 
ter validity of a mean-field treatment. 



IV. CONCLUSION 

In conclusion, we have presented an effective theory 
for the merging of two vortices advected by the effective 
noise of the other vortices. Several quantities including 
the merging time r„j and the pdf P{r, t) have been de- 
rived analytically and are in excellent agreement with nu- 
merical simulations. The results of this study point to the 
weak efficiency of two-body processes, as the fast rota- 
tion of the vortex pair averages out the effective advecting 
noise due to the other vortices. We numerically find that 
a similar effective theory, where three vortices of different 
sign now participate to the merging process, leads to a 



dynamics would still be treated exactly, but in the effec- 
tive noisy background generated by all the other vortices. 



Appendix A: Conservation of the vorticity peak 

It is numerically observed that the vorticity peak is al- 
most conserved in 2D decaying turbulence [l^ , an impor- 
tant ingredient of the scaling theory. In this Appendix, 
we would like to motivate this claim in terms of statistical 
mechanics and kinetic theory of the 2D Euler equation. 

The Miller-Robert-Sommeria (MRS) statistical theory 
[36l . [37| is relatively well-suited to describe the rapid 
merging of two vortices. Knowing the initial condition 
and assuming crgodicity (efficient mixing) the statisti- 
cal theory predicts the shape of the new vortex result- 
ing from this merging. Furthermore, the mixing process 
can be described in terms of relaxation equations [ssj 
obtained from a maximum entropy production principle 
(MEPP). Assuming, for simplicity, that the initial vor- 
tices have a uniform vorticity w = ctq, the relaxation 
equations can be written 



did 

'dt 



+ u • Vcj = V • [^(r, t) (Vw + l3{tp{aQ - oJ)\/i')] 



m 



J D{r,t)Vuj -Vtjjdr 
' J D{v, t)uj{ao - w)(VV')2 dv ' 



(Al) 



(A2) 



where cJ(r, t) is the coarse-grained vorticity field and 
D{r,t) the "turbulent" diffusion coefficient. If D{r,t) is 
constant, the relaxation equations converge towards the 
statistical equilibrium state given by a Fermi-Dirac-like 
distribution 



ZU(r) 



1 _|_ g<To(^*iA(r)+Q) ■ 



(A3) 



However it has been observed at several occasions (see, 
e.g., [l^liOl) that relaxation is incomplete due to lack of 
ergodicity. This is clearly a limitation of the statistical 
theory. In order to take incomplete relaxation into ac- 
count, it has been proposed that the diffusion coefficient 
_P(r, t) varies in space and time. Qualitative arguments 
[4l| . or more sophisticated kinetic theory [42|, indicate 
that the diffusion coefficient should be proportional to 



faster decay of the density, n{t) t \ after a transient ^-^c local variance uj2 = lu^-lu^ = lu{(To - l^) of the vor- 



regime where n(t) ~ We therefore conclude that 

three-body processes should dominate, at least at small 
vortex density. This result can be understood within a 
simple mean-field kinetic theory describing the merging 
process as the collision of a ballistic dipole with an iso- 
lated vortex. This same theory applied to geostrophic 
turbulence in three dimensions leads to n{t) ~ t'^^'^, in 
good agreement with numerical simulations. It would 
be certainly interesting to generalize this study for an 
increasing number iV > 4 of tagged vortices, yet small 
enough for the numerical simulations to reach very small 
effective densities {i.e. small Tc). The N vortex mutual 



ticity ffuctuations. It must also decrease in time as the 
fluctuations decay. Therefore 



D{r,t) ^ A{t)iJ{r,t)[ao -uJ{r,t)], 



(A4) 



where A{t) tends to zero on a few dynamical times. In the 
"mixing zone" where the fluctuations are strong, the dif- 
fusion coefficient is large and the system rapidly reaches a 
Fermi-Dirac-like distribution (|A3[) . with an inverse tem- 
perature /?' and a chemical potential a' that may be dif- 
ferent from the equilibrium ones. On the other hand, 
from Eq. (|A4I) . one sees that the diffusion coefficient 



7 



is small in the regions where the fluctuations are weak. 
This concerns the core of the vortex (aJ(r, t) ao) where 
the vorticity is maximum and the region surrounding the 
vortex {Lo{r, t) 0) which is not sampled by vortic- 
ity. Therefore, the relaxation towards the Fermi-Dirac- 
like distribution (jA3p will be very slow in these regions. 
Since, in parallel, the diffusion coefficient decreases with 
time (because of the term A{t)), the system will not have 
time to reach the Fermi-Dirac-like distribution in these 
regions. Therefore, mixing will be relatively inefficient in 
the core and at the periphery of the vortex. These ar- 
guments of kinetic theory may explain why the vorticity 
peak is relatively well conserved during a merging and 
why the vortex is more confined than what the statis- 
tical theory predicts assuming complete mixing. These 
arguments have been developed in Sec. 4.3 of |43|. 



ordinates. Equations (|9I10|) are therefore rewritten as 
dx y 



dy X 

- = ^-yV. + xv,. 



(Bl) 



(B2) 



Introducing the convenient notations r = (x, y) and w = 
('7q,?7^), the Fokker-Planck equation for P(r, w,t) is 



dP 
'dt 



d_ 

dx 



d \ ( X 
dy 



y 

P + xj]a + yvp 1 P 
a f dp 

dw \ dw 



Pw .(B3) 



The stationary solution {dtP — 0) of this Fokker-Planck 
equation is obtained when the r.h.s. and the l.h.s. vanish 
individually. The cancelation of the r.h.s. implies that 



Ps{r,Mv) ^ C(r)e 



(B4) 



Inserting this expression in the l.h.s. of the Fokker- 
Planck equation, we obtain C(r) = constant. Therefore, 



Appendix B: The Fokker-Planck equation 



In order to write down the Fokker-Planck equa- 
tion associated with the stochastic process defined by 
Eqs. (|9ll0ll3p . it may be simpler to use cartesian co- 



Ps(r, w) cx e' 



(B5) 



It is noteworthy that the steady state of this Fokker- 
Planck equation is spatially homogeneous, i.e. Ps{j^) = 
constant, and Psir) = r in radial coordinates. 
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